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By combining a recent estimate of the total B — » X s "f branching fraction at 0{a 2 a ) with a detailed 
analysis of the effects of a cut _E 7 > 1.6 GeV on photon energy, a prediction for the partial B — > X s -y 
branching fraction at next-to-next-to-leading order in renormalization-group improved perturbation 
theory is obtained, in which contributions from all relevant scales are properly factorized. The result 
Br(_B — > X s j) = (2.98 ± 0.26) • 10~ 4 is about 1.4c lower than the experimental world average. This 
opens a window for significant New Physics contributions in rare radiative B decays. 



I. INTRODUCTION 

The inclusive decay B — > X s ^ is an important exam- 
ple of a flavor-changing neutral current process, which 
has been used to test the flavor sector of the Standard 
Model. Many groups have worked on improving the the- 
oretical analysis of this process so as to keep pace with 
refinements in the measurements of its branching frac- 
tion. The effective weak Hamiltonian at next-to-next-to- 
lcading order (NNLO) has been obtained by calculating 
multi-loop matching coefficients and anomalous dimen- 
sions |H [J [1 While the fermionic NNLO corrections 
to the b — > 57 matrix elements have been known for some 
time 0, complete NNLO corrections are presently only 
available for the electro-magnetic dipole operator [a, [3 ■ 
However, an approximate result for the NNLO charm- 
penguin contributions has just been published Q. Com- 
bining these ingredients, a first estimate of the B — * X s j 
branching ratio at NNLO has been presented in 

A complication in the analysis arises from the fact 
that measurements of the B — > X s ^ branching fraction 
impose stringent cuts on photon energy (defined in the 
B-meson rest frame), E~ t > Eo, with Eq in the range 
between 1.8 to 2.0 GeV. The standard treatment is to 
extrapolate different measurements to a common refer- 
ence point Eo = 1.6 GeV using phenomenological mod- 
els LLfJ. In that way, the experimental world average 
Bi(B ^_^7) = (3.55 ± 0.24™ ±0.03) -l(r 4 has been 
derived nil. The first error is statistical, the second one 
systematical, the third one is due to the extrapolation 
from high Eq to the reference value, and the last error 
accounts for the subtraction of B — > Xd"f background. A 
theoretical result for the branching ratio with a cut at 
Eq = 1.6 GeV has been derived in [!j using two-loop cal- 
culations of the photon-energy spectrum in fixed-order 
perturbation theory [l2l ll^ . It has been argued that 
the extrapolation from the total to the partial branching 
fraction does not introduce additional theoretical uncer- 
tainties. This assertion is questionable because of the dy- 
namical relevance of a soft scale A = m^—2EQ sa 1.4 GeV, 
whose value is significantly lower than the 6-quark mass. 

Accounting for the photon-energy cut properly re- 
quires to disentangle contributions associated with the 



hard scale fih ~ fnbi the soft scale /io 
mediate scale ^ 



A, and an inter- 



\fmbK set by the typical final-state 
hadronic invariant mass. When the cut value Eq is cho- 
sen sufficiently low, A becomes a short-distance scale, 
and renormalization-group (RG) improved perturbation 
theory can be employed to calculate the effects of the 
photon-energy cut using a multi-scale operator product 
expansion [l4[ . We have recently performed a systematic 
analysis of these effects at NNLO. Two- loop corrections 
at the soft scale were calculated in |l5J, while those at 
the intermediate scale were computed in |l6j . Here, the 
analysis is completed by extracting the two-loop hard 
matching corrections from a comparison with fixed-order 
calculations of the photon spectrum. 

Using this method, we compute the fraction of all 
B — > X s j events with E 1 > 1.6 GeV with a perturbative 
precision of 5%. At this level of accuracy several other, 
nonperturbative effects need to be evaluated carefully. 
The event fraction receives hadronic power corrections 
~ (Aqcd/A)" governed by .B-meson matrix elements of 
local operators. The leading correction (n = 2) is known 
and turns out to be small, but terms with n > 3 are 
presently unknown. Recently, a new class of enhanced 
A-QCD/mb corrections to the B — > X s ^ decay rate has 
been identified, which involve matrix elements of nonlo- 
cal operators |l7j| . A model analysis using the vacuum 
insertion approximation indicates that these corrections 
affect the total decay rate at the level of a few percent. 

Combining our result for the event fraction with the 
prediction for the total branching fraction from [8J, |9| , we 
obtain 



Br(B -> X sl ) = (2.98 ± 0.26) • 10^ 



(1) 



for Eq = 1.6 GeV, where we have added in quadrature 
the uncertainties from higher-order perturbative effects 
(^4%), hadronic power corrections (5%), parametric de- 
pendencies (4%), and the interpolation in the charm- 
quark mass (3%). Two-loop perturbative corrections at 
the intermediate and soft scales significantly lower the 
branching fraction with regard to the fixed-order result 
given in |8J| , and they increase the theoretical uncertainty. 
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II. SCALE SEPARATION AND RESUMMATION 

At leading power in Aqcd/w(,, the B — * X s j decay 
rate with a cut on photon energy obeys the factorization 
formula 



T{Eo) 
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where p + = m b 
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mi, — 2Eq. The function 



contains hard quantum corrections, the jet function 
J describes the physics of the hadronic final-state jet, 
and the shape function S parameterizes bound-state ef- 
fects inside the S-meson In the region of interest 
to our analysis, in which the quantity A sa 1.4 GeV can 
be treated as a perturbative scale, the double convolu- 
tion inte gral in Q can be evaluated using short-distance 
methods [lj] . In order to isolate the effect of the photon 
cut, we focus on the fraction of all B — > X s j events that 
pass the cut E y > E , defined as F(E ) = T(E )/T(0). 
Since the total decay rate can be computed in fixed-order 
perturbation theory at the hard scale fih, the event frac- 
tion obeys a factorization formula of the same form as J2J), 
but with a different hard function h in place of |i? 7 | 2 . 

Contributions associated with different mass scales can 
be separated from each other by evolving the various ob- 
jects in the factorization formula away from a common 
renormalization scale /i to different "matching scales" , 
where they can be calculated reliably using fixed-order 
perturbation theory. In this process, single and double 
logarithms of ratios of the different scales are resummed 
to all orders. The matching scales should be taken close 
to the default values fih = m b for _ff 7 , fii — y'm&A for J, 
and /io = A for S. An elegant expression describing the 
RG evolution of the jet function was derived in [2(|. It 
involves an associated jet function j, which is related to 
J by a Laplace transform. The same technique can be 
applied to describe the evolution of the shape function 
in terms of an associated soft function s. Inserting these 
results into @, the integrations over p + and w can be 
performed, leading to 

F(E ) = U{nh,m,Ho\y) 
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+ SF(E ), (3) 



where rj = 2av{Hi, Mo) > 0. Despite appearance, this re- 
sult is independent of the choice of /i and of the three 
matching scales fi^, /i,, and ^o- F° r the special case 
fi = Hi, relation coincides with a formula derived pre- 
viously in pH . The terms in the first two lines arise from 



the RG resummation of single and double logarithms. 
The precise form of the evolution factor U will be given 
in |22j. The exponent 
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is an integral over the cusp anomalous dimension. The 
functions h, j, and ^in the third line contain the match- 
ing corrections at the hard, intermediate, and soft scales, 
£espectively. The two-loop expression for the jet function 
j has been obtained in , while the two-loop result for 
the soft function s can be deduced from [l^ . In the ar- 
gument of these functions, d n means a derivative with 
respect to the quantity 77. The hard function can be de- 
rived by matching J3J with the fixed-order expression for 
the photon spectrum derived in [T2I Il3| . At the default 
scale /j,h = ra b , we find 
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where the term with i = j = 7 is to be excluded from 
the sum. The one-loop result agrees with ^4[, and the 
two-loop coefficients are given by 
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(6) 



The constants S& w 1.216 and Sn a « —4.795 have been 
obtained by numerical integration in [l2j . The complete 
expression for h(m{,/ fih) including scale dependence will 
be given in |22l |. 

In the last line of (J3J, the hadronic quantity /z 2 pa- 
rameterizes the -B-meson matrix element of the kinetic 
operator in heavy-quark effective theory. The ellipses 
represent unknown hadronic power corrections of order 
(Aqcd/A) 3 . The polynomial 



P 3 (S) 
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(7) 



derives from the (mj, — p + ) 3 prefactor in (0. The term 
SF(Eq) in © collects terms that are power suppressed in 
the ratio 5 = A/m& ~ 0.3. These corrections are known 
in fixed-order perturbation theory, but no resummation 
is available. We express them as a perturbative series in 
powers of a s (/x). The one-loop result for SF has been 
given in |T^. while the dominant two- loop contribution 
has recently been calculated in [T2I IT3I ] . For subleading 
contributions from other operators in the effective Hamil- 
tonian only the fermionic two-loop corrections are known 
|23| : they have a negligible impact on our results. 
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FIG. 1: Scale dependence of the ratio T. The curves show the 
effect of varying the hard scale fih (dotted), intermediate scale 
/Xi (dashed), soft scale fio (dash-dotted), and the reference 
scale jj, (solid) by a factor of V2 about their default values. 
In the left plot the resummation is performed for the leading- 
power terms only, while in the right plot the P3 term is also 
included (see text for further explanation). 



III. RESULTS AND CONCLUSIONS 

In ©, the parameters mt and n\ are defined in the 
on-shell scheme. To improve the perturbative behav- 
ior, one should eliminate them in favor of appropriately 
defined short-distance quantities. We use the "shape- 
function scheme" proposed in [24). in which low-scale 
subtracted heavy-quark parameters are defined via the 
moments of the renormalized i?-meson shape function, 
regularized with a hard cutoff /i/. The two- loop rela- 
tions between the shape-function and pole-scheme pa- 
rameters have been derived in [^J. We adopt the con- 
ventional choice fif = 1.5 GeV and use the values nib = 
(4.61 ± 0.06) GeV and n\ = (0.15 ±0.07) GeV 2 extracted 
from a fit to moments of inclusive .B-decay spectra [2l| . 

Ref. does not provide a result for the total B —> X s j 
branching fraction. The most inclusive quantity consid- 
ered is the partial branching fraction defined with a mild 
cut at Eq — 1 GeV. In order to combine this prediction 
with our RG-improved result for the event fraction, we 
define the ratio T = F (1.6 GeV) /F (1.0 GeV). We eval- 
uate the various matching scales at their default values, 
namely = mj, fj,o = A, and fa = fi = ^/m^A, where 
the values of A are different in the numerator and de- 
nominator. In order to study the residual scale depen- 
dence, we then vary each of the four scales by a factor 
of \/2, correlated between numerator and denominator. 
While the ratio T is formally independent of these scales, 
the residual dependence of the truncated perturbative 
expression can be taken as an estimate of higher-order 
effects. The results are depicted in Figure^ Not surpris- 
ingly, the dominant effect arises from varying the lowest 
scale, fiQ r-j A, while the scale variations at the interme- 
diate and high scales have a lesser impact. 

The term proportional to ^(A/m?,) in @ includes a 
subset of power corrections associated with a phase-space 
factor. While it was possible to perform the scale separa- 



tion for these terms, treating them in a different way than 
the remaining power corrections in SF(Eq) is somewhat 
arbitrary. The right plot in the figure refers to the form of 
the factorization formula shown in J3J), while the left plot 
corresponds to expanding out (p^ — 1) in fixed-order per- 
turbation theory and including it in the SF term. While 
the stability with respect to variations of the soft scale 
is better in this case, the perturbative corrections turn 
out to be smaller when performing the resummation for 
the P3 term. The shift in central value between the two 
schemes is about 3%, which is inside the error bar. This 
effect hints at the importance of RG resummation for the 
power corrections. 

To quote our final result we take the average of the 
two schemes and assign an asymmetric error reflecting 
the scale variation. This yields 



T = 0.93l°;°5Pert ± 0.02 hadr ± 0.02 pa 



(8) 



The event fraction F(Eq) receives hadronic power cor- 
rections not suppressed by inverse powers of mj, but 
only by powers of the soft scale A [14). These correc- 
tions are governed by .B-meson matrix elements of lo- 
cal operators. The leading effect proportional to /i 2 in 
© is small mainly due to the smallness of its coefficient 
77(1 — r?)/6. Generically, we expect subleading correc- 
tions to scale like ^(Aqcd/A) 3 , for which we assign a 
2% uncertainty. The main uncertainties from parame- 
ter variations are tO.9% for a s (m z ) = 0.1189 ± 0.0020, 
±0.4% for m c /m b = 0.26 ± 0.03, and ±0.1% for m b = 
(4.61 ± 0.06) GeV. We also include a variation of ±1.2% 
due to the fact that the three-loop anomalous dimension 
of the shape function is yet unknown p2| . Our value for 
the ratio T is lower than the estimate T = 0.963 obtained 
in fixed-order perturbation theory [9( . Moreover, we find 
that there is a significant theoretical uncertainty inherent 
in the calculation of T. 

In order to complete the analysis we need as input 
the theoretical result for the B — ► X s j branching frac- 
tion with Eq = 1.0 GeV, which we take from the fixed- 
order NNLO calculation of 0, These authors find 
Br(_B -> X sl ) = (3.27 ± 0.23) • 10" 4 for E = 1 GeV, 
where the error has been obtained by adding in quadra- 
ture the uncertainties from higher-order perturbative cor- 
rections (3%), nonperturbative effects (5%), parameter 
dependencies (3%), and the interpolation in the charm- 
quark mass employed in the NNLO estimate of charm- 
penguin loop graphs (3%). The B — > X s j branching 
fraction receives incalculable power corrections starting 
at order Aqcd/ 'tij, which cannot be described using the 
operator product expansion 26] . While we disagree with 
|8| on the statement that these effects are suppressed by 
a power of a s (mt>), we nevertheless believe that the cor- 
responding uncertainty should not be larger than 5%. 
Recently, a new class of enhanced power corrections has 
been identified ^3 ■ At tree level their effects are parame- 
terized in terms of _B-meson matrix elements of nonlocal 
four-quark operators. Using the vacuum insertion ap- 
proximation, a reduction of the total branching fraction 
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between 0.3% and 3% has been found. Accounting for 
this effect lowers the central value from 3.27 to 3.22. 

Our final prediction for the B — > X s j branching frac- 
tion with Eq = 1.6 GeV is 



Br(B -> X sl ) 

= C2.98±g;i| pert ± 0.16 hadr ± 0.11 pars ± 0.09, 



(9) 

•io- 4 , 



where we have combined errors of the same type in 
quadrature. This appears justified, since theoretical cor- 
relations in the calculations of the total branching frac- 
tion and the event fraction F(Eo) are small. In Q we 
have combined all uncertainties in quadrature. A more 
conservative approach would be to add the errors lin- 
early, in which case the error becomes • Compared 
with the result Br(B -> X s j) = (3.15 ± 0.23) • 10~ 4 ob- 
tained in our central value in © is lower by about 
5% and, more importantly, the perturbative uncertainty 
is larger by almost a factor of 2. Both changes are a re- 
sult of significant two-loop corrections encountered at the 
intermediate and soft scales, ^ ~ and /iq ~ A. 

Our theoretical prediction for the B — ► X s j decay rate 
is consistent with the present experimental world aver- 
age, as is reflected in the ratio 



Br{B - X sl ) 



cxp 



Bt(B -> X sl ) th 



= 1.19 ±0.09 cxp ± 0.10th 



(10) 



However, whereas for a long time the experimental re- 
sult used to be lower than the theoretical one, it is now 
about 1.4 standard deviations larger. Since in many ex- 
tensions of the Standard Model the contributions from 



New Physics are expected to interfere constructively with 
the Standard Model b — ► S7 amplitude, the situation 
has changed from one where New Physics models were 
rather tightly constrained to one where there is now room 
for speculation about how the central number in 110|) 
could be explained in terms of loop contributions con- 
taining new heavy particles. Consider, e.g., the case of 
type-II two-Higgs-doublet models. Whereas in the past 
there used to be bounds on the charged-Higgs mass from 
B — > X s j on the order of 500 GeV, the ratio JTJJ could 
now be explained with a New Physics contribution from 
a charged Higgs in just that mass range. 

In summary, we have presented the first NNLO pre- 
diction for the B — ► X s j branching fraction in which 
the effects of a photon-energy cut E~ ( > 1.6 GeV have 
been properly taken into account. Low-scale perturba- 
tive corrections lower the prediction for the branching 
ratio and introduce a significant theoretical uncertainty 
even at NNLO. Our result is about 1.4 standard devi- 
ations lower than the world-average experimental value. 
This re-opens the door for explorations of New Physics 
contributions to rare flavor-changing _B-decay processes. 
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